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A body- contouring method for alleviating at subsonic speeds the 
adverse interference at the root of a high-aspect-ratio sweptbaek wing 
was investigated. Several bodies used in combination with a 35 ° swept- 
back wing were modified, as proposed by Kuchemann, by shaping the body 
so as to counteract the distorting velocities at the wing-body junction. 

Although modifying the body shape did not significantly affect the 
aerodynamic characteristics at subcritical speeds, beneficial results 
were obtained at free-stream Mach numbers above the critical. Improved 
aerodynamic characteristics were evidenced by large reductions of drag, 
an increase in lift- curve slope, and a reduced change of pitching-moment- 
curve slope with increasing Mach number. 

Additional tests, which involved changes in wing position relative 
to the body, indicated that wing position had little effect on the aero- 
dynamic characteristics of the unmodified wing-body combinations. 


INTRODUCTION 


The effect of sweep for delaying the onset of compressibility drag 
has generally been somewhat less beneficial than indicated by simple- 
sweep theory. This is caused, at least in part, by an adverse velocity 
distribution at the root of the swept wing. The subsonic wing theories 
of references 1 and 2 indicate that the velocity distribution at the root 
of a swept wing is, because of reflection at the plane of symmetr y, dis- 
torted in such a manner that the chordwise pol 
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is displaced rearward. Furthermore, the maximum velocity at the root of 
the swept wing may he considerably higher than that for an infinite yawed 
wing and, consequently, the premature formation of strong local shocks 
in the vicinity of the root can be expected. 

The velocity distortion can be alleviated by altering either the 
wing geometry near the root or by contouring the body in the vicinity of 
the wing-body junction. Shaping the sides of the body to conform to the 
general shape of a streamline on a yawed wing has been suggested in 
reference 3* Experimental investigations described in references 4 and 5 
were based on this method and some beneficial results were obtained. A 
modified swept-wing-root profile was also investigated in the research 
reported in reference d. 

The design of wing-fuselage junctions for subsonic speeds has been 
given an exhaustive treatment by Kuchemann (refs. 6 and 7) • In one of 
the methods proposed by Kuchemann, the body is represented by a cylinder 
on which ring vortices are distributed so that the induced axial compo- 
nent of velocity cancels the swept-wing interference velocity. By inte- 
gration of the induced lateral velocity, the radial modification neces- 
sary to shape the wing-body junction is determined. An experimental 
investigation at low speeds (ref. 8) has indicated the efficacy of this 
approach. 

The present report gives the results of an experimental investiga- 
tion of this method at high subsonic speeds. A summary of the theoreti- 
cal method for determining the required body modifications and a detailed 
description of the application of this method as used in the present 
investigation are presented in the Appendix of this report. Although the 
present investigation was conducted primarily to evaluate a particular 
body contouring method, several related studies were conducted. In order 
to obtain information concerning the effect of the swept-wing position 
relative to the body on the aerodynamic characteristics, a body was 
tested with the wing in forward and in rearward positions. A body was 
also tested with the wing at the center line of the body and in a raised 
position. Additional information concerning the effect of mutual inter- 
ference between body- and wing-induced velocity fields was obtained by 
testing the wing in combination with a body having a cylindrical mid- 
section. 

It should be recognized that the present method of body contouring 
is primarily concerned with eliminating the interference at the root of 
a swept wing. The present method should not be confused with the methods 
of. references 9 a nd 10, which are based on altering the axial distribu- 
tion of cross-sectional area to minimize the wave drag at Mach numbers 
of 1 or greater. 


iONFIDENTIAS* 



NACA RM A54A22 


CONFIDENTIAL 


NOTATION 


A 

b 

C D 

c L 

Cm 


aspect ratio 

wing span 

drag coefficient, 

lift coefficient. 


drag 

“is" 

lift 

qS 


pitching-moment coefficient, 
local wing chord 

wing mean aerodynamic chord. 


pitching moment about c/4 


qSc 


f b / 2 „2 


c dy 


jb / 2 


c dy 


Cj chord at wing-body junction (chord through the point of intersec- 
tion of basic body and midchord line of the swept wing) 

d G maximum body diameter 

Z body length (distance from nose to theoretical point of closure) 

M free -stream Mach number 

M cr wing critical Mach number based on simple -sweep concepts 
Mdes design Mach number 
P pr e s sure c oef f i c i ent y 

(local static pressure) - (free-stream static pressure) 

- 

1 2 

q free-stream dynamic pressure, — PV 0 

r body radius 

r Q maximum body radius 

r j radius of basic body at intersection with midchord line of swept 
wing 
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ving area 

wing thickness- to-chord ratio 

streamvise perturbation velocity- 
free- stream velocity 
distance behind body nose 
lateral distance from model center line 

ordinate of wing surface, dimensionless with respect to the local 
wing chord 

vertical distance above model center line 
angle of attack 
Vl - M 2 

^ 1 ' M des 2 . 
wing taper ratio 

angle of sweep, positive when swept back 

effective angle of sweep, tan A e = ^ an — 

^des 

distance behind the leading edge of the wing -body .junction, dimen- 
sionless with respect to the wing chord at the wing-body junction 

free-stream mass density 

radial modification associated with the vortex cylinder 
modification to basic body radius in horizontal plane of symmetry 
modification to basic body radius in vertical plane of symmetry 
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APPARATUS AND MODELS 
Apparatus 


The tests were conducted in the Ames l6-foot high-speed wind tunnel 
with the models mounted on a sting support as shown in figure 1. The 
normal and chord forces and the pitching moment were measured with elec- 
trical strain gages enclosed within the model. Multiple -tube mercury- 
manometers, connected to pressure orifices in the model by flexible 
tubing, were photographed to provide records of the pressure distribution 
on the model surface. 


Models 


In this report, for reasons of brevity, those models designed with- 
out consideration of the wing-body interference problem are designated 
as "basic" and those designed with consideration of the interference 
problem are designated as "modified." The wing used in combination with 
the various bodies had an aspect ratio of 6, a taper ratio of 0 . 5, and 
NACA 6^A015 sections normal to the 50-percent-chord line, which was 
swept back 35° • (See fig. 2.) For all the wing-body combinations tested, 
the center lines of the bodies were located in the chord plane of the 
wing, except for the model shown in figure 2(c) which was tested with the 
wing chord plane at the body center line and also with it 16.1 percent of 
the maximum body diameter above the body center line. 


Configurations having basic bodies .- The bodies of revolution shown 
in figures 2(a), 2(b), 2(c), and 2(d) were shaped in accordance with the 
Sears -Haack formula 


r 


r o 



except for the afterportion of the model shown in figure 2(b) where 
straight-line elements were used as shown. 


The model shown in figure 2(e) had a body with a cylindrical mid- 
section approximately twice the length of the wing root chord. The body 
radii forward and aft of the cylindrical midsection were computed using 
the Sears -Haack formula for which the corresponding values for z/2 were 
taken as the body length forward and aft of the cylindrical midsection, 
respectively, with r 0 equal to the cylinder radius. 


All the bodies were truncated, as indicated in figure 2, in order 
to permit mounting on the sting. The fineness ratios of the bodies were 
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computed using the theoretical length of the body to closure and the 
maximum body diameter. 

Configurations having modified bodies .- The bodies chosen for modi- 
fication were those shown in figures 2(a), 2(c), and 2(d). These bodies 
were contoured in the vicinity of the wing-body junction so as to mini- 
mize the interference velocities at the design Mach number, 0.87> which 
was equal to the predicted critical Mach number of the swept wing (based 
on simple-sweep concepts). 

A summary of the method for calculating the body modifications is 
given in the Appendix. The resultant bodies were not bodies of revolu- 
tion in the vicinity of the wing-body junction but were elliptical in 
cross section. Details of the contouring are presented in figure 3* 

(See also the photographs in fig. 1.) 

The basic bodies chosen for modification were also tested without 
the wing in order that the predicted body pressures, used in the contour- 
ing calculations, could be compared with measured body pressures. 


TESTS AND PROCEDURE 


The models were investigated through a Mach number range from 0.50 
to 0.9b, with a corresponding Reynolds number variation (based on the 
wing mean aerodynamic chord) from 2.6xl0 6 to 3*^Xl0 6 (fig. 4) . The test 
data have been corrected for the effects of blockage and tunnel-wall 
constraint by the methods of references 11 and 12. 

The base drag, computed from the base areas of the bodies and the 
difference between measured base pressures and the free-stream static 
pressure, has been subtracted from the drag measurements. The drag data 
are believed to be slightly in error due to an interaction of the normal 
force and pitching moment on the chord-force measurements of the balance. 
This interaction was known to be small and inconsistent so no correction 
could be applied. 


RESULTS AND DISCUSSION 
Force Studies 


The lift, drag, and pitching-moment data are shown in figure 5 for 
the nine configurations. Cross plots of these data, which summarize the 
drag, lift, and static longitudinal stability of the models, are pre- 
sented in figures 6 to 9- 
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Drag .- It is apparent from the basic data shown in figure 5 that 
the drag for the configurations having modified bodies was much less 
than that for the basic configurations at Mach numbers above the criti- 
cal (the predicted critical Mach number for the swept wing, based on 
simple-sweep concepts, is approximately 0.87) . Modifying the bodies 
also resulted in reductions in drag at the higher lift coefficients for 
free-stream Mach numbers slightly less than O. 87 . The zero-lift drag 
at subcritical Mach numbers was not greatly affected by the contouring. 

The body modification increased the M^ch number for drag diver- 
gence, defined as the Mach number for which 



by approximately 0.02. (See fig. 6.) More striking, however, was the 
reduction in drag above the divergence Mach number, the reduction in 
zero-lift drag at a Mach number of 0.9 1 *- being of the order of 40 percent. 

The drag for the wing in the rearward position was only slightly 
less than that for the wing in the forward position (fig. 6(a)). The 
configuration having a body with cylindrical midsection, when compared 
with the one using the basic body of the same size, had less drag at 
Mach numbers above the critical but slightly greater drag at subcritical 
speeds (fig. 6(b)). Raising the wing did not affect the zero-lift drag 
(fig. 6(c)). The same effects are shown for a lift coefficient of 0.3 
in figure 7. 

Lift .- The variation of lift-curve slope (evaluated at zero lift) 
with Mach number is shown in figure 8. Although minor differences in 
lift-curve slope for the various basic and modified models are evident 
at low Mach numbers, the most significant effects occurred at Mach num- 
bers near and above the critical. Modifying the body shape resulted in 
less variation of lift-curve slope with Mach number and an increase in 
the lift-curve slope at supercritical Mach numbers. 

A comparison of the data for the models with the wing forward and 
with the wing aft (fig. 8(a)) indicates that the lift-curve slope for 
the model with the wing aft rose to a somewhat lower maximum value near 
the critical Mach number but then decreased in a somewhat similar manner 
with increasing supercritical Mach number. The model having a body with 
a cylindrical midsection had a slightly lower lift-curve slope at low 
speeds than the basic body of the same size (fig. 2(d)). The maximum 
slope was reached at a Mach number slightly greater than that for the 
basic body, but the rate of change of lift-curve slope with increasing 
supercritical Mach number was greater (fig. 8(b)). Raising the wing to 
a position above the center line of the body did not affect the lift- 
curve slope (fig. 8(c)). 
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Moment .- The variation of the static longitudinal stability 

f dc m \ 

parameter, I - ) , with Mach number is presented in figure 9« 

° L c L =o 

Modifying the body shape resulted in smaller changes in the stability 
parameter with changes in Mach -number. 

The longitudinal position of the wing with respect to the body had 
a large effect on the stability parameter (fig. 9( a ) ) • The change in 
stability occurring with increasing Mach number was most pronounced with 
the wing in a rearward position. The effect of replacing the midsection 
of the body by a cylinder was relatively unimportant (fig. 9(b) )• Rais- 
ing the wing resulted in small changes in the stability parameter 
(fig. 9(e)). 


Pressure Studies 


The effects of the body modification on the pressure distribution 
over the inner half of the wing at zero angle of attack are illustrated 
in figure 10 for three high subsonic Mach numbers. With the basic bodies 
the wing pressures near the root were less negative than those along the 
outer panel of the wing, but, as predicted by theory, an unfavorable 
distortion of the pressure distribution occurred near the root, the loca- 
tion of the point of minimum pressure being shifted rearward. 

The bodies had been modified in an attempt to obtain the same pres- 
sure distribution near the root chord as over the outer portion of the 
wing. Although this was far from being accomplished, a more favorable 
distribution of pressure over the root was obtained. Specifically, the 
point of most negative pressure was moved from a position behind the 
root midchord to a position forward of the midchord. The influence of 
the body modifications extended to a considerable distance from the root, 
and the absence of severe pressure gradients near the trailing edge of 
the inner portion of the wing is particularly noteworthy. 


CONCLUDING REMARKS 


These tests show that a considerable improvement in the aerodynamic 
characteristics of wing-body combinations, employing a swept wing of 
large thickness- to- chord ratio and large aspect ratio, can be obtained 
at high-subsonic Mach numbers by contouring the body in the vicinity of 
the wing-body junction according to the ring-vortex design method sug- 
gested by Kiichemann. This type of body modification did not significantly 
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affect the aerodynamic characteristics at Mach numbers below the critical 
of the swept wing. At Mach numbers above the critical, improved aero- 
dynamic characteristics were evidenced by large reductions of drag, an 
increase in lift-curve slope, and a reduced change of pitching-moment- 
curve slope with increasing Mach number. Vertical and longitudinal 
changes of the position of the wing on the body had little effect on the 
aerodynamic characteristics of the unmodified wing-body combinations. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Jan. 22, 195^ 
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APPENDIX 

METHOD OF BODY MODIFICATION 


The design objective is to alter the body shape so that the velocity 
distribution at zero lift, for a given subsonic Mach number, will be the 
same at the wing-body junction as that on an infinite yawed wing at the 
same Mach number. Computations of the required body contours are based 
on a method proposed by Kuchemann (refs. 6 and 7) in which linearized 
theory is used throughout. 


Perturbation Velocities at the Wing-Body Junction 


The streamwise perturbation velocity uj . in the plane of the wing 
and at the wing-body junction is regarded as being composed of the fol- 
lowing additive components: 





u c u m 



(Al) 


where 

Ub perturbation velocity induced by a body alone 

u w perturbation velocity of an infinite yawed wing 

u c distortion velocity at the root of a swept wing 1 

u m perturbation velocity induced by modifications of the basic body 
shape 

In order that the velocity at the wing-body junction be the same as 
for an infinite yawed wing, it is necessary that the body modification 
result in a perturbation velocity of sufficient strength to cancel the 
sum of the body perturbation velocity u^ and the wing-root distortion 
velocity u c } that is, the condition 



is to be satisfied along the wing-body junction (reference chord). 
^-In this procedure the body is considered to be a reflection plane. 
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Design Mach Number 


In general, the compressibility corrections will be different for 
the two velocity ratios u c /V 0 and u^/v 0 to be canceled by means of the 
body modification. Because the variation with Mach number is different 
for these two velocity ratios, it is necessary to select a design Mach 
number for use in computing the body modification. The most beneficial 
results of the body contouring are expected at supercritical speeds; 
therefore, it seems rational to use the critical Mach number of the yawed 
wing as the design Mach number, since this is the highest subsonic Mach 
number for which linearized theory can be expected to give reasonable 
results. The NACA 6^015 airfoil used has a theoretical critical Mach 
number of O .87 when yawed 35°. 


Body-Induced Velocities 


The body-induced velocities may be computed by means of any of the 
methods proposed for bodies of revolution in axially symmetric flow. The 
slender-body approximation of Laitone (ref. 13) is used here. 2 The for- 
mula is 



(n) 

where S(x) is the cross-sectional area of the body and S v (x) is the 
nth derivative of this area with respect to x. The body equation for 
use in equation (A3) is, for the bodies of this report, the Sears-Haack 
formula 


(A3) 


r 



K i -?)T 


(A4) 


It should be noted that equation (A3) includes the effect of com- 
pressibility. The first two terms are sufficient to give a reasonably 

2 The use of Laitone' s method, which depends on expanding the expression 
for the cross-sectional area of the body in a power series, is permis- 
sible here because the bodies used are described by an analytically 
continuous function. 
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good prediction of the body-induced velocities near the center section 
of slender bodies. A comparison of calculated and measured body pres- 
sures is presented in figure 11. 


Distortion Velocity at the Root of a Swept Wing 


The velocity near the root of a swept wing is known to be distorted 
from that of an infinite yawed wing. Klichemann has determined the magni- 
tude of the distorting velocity at the root of an untapered wing of 
infinite aspect ratio in incompressible flow to be 


a az(S) 

cos A 

d£ 1 + sin A 

= in 

Vo n 1 - sin A 


(A5) 


where Z(l) is the ordinate of the airfoil surface, made dimensionless 
with respect to the wing chord, and A is the angle of sweep (positive 
when swept back) . 

For compressible flow, the three-dimensional form of the Prandtl- 
Glauert transformation (ref. l4) is applied to equation (A5) . At a 
given design Mach number, M^ es , equation (A5) becomes 


. . cos A e — 

/ He ( £ ) \ _ 1 + sin Ae 

V Vo / de s *3des 1 " sin A e 


(A6) 


where 

%es ~ ^ ^des 


and the effective angle of sweep is related to the actual sweep angle 

by 


, A tan A 
tan Aq = 

%es 


The distortion velocity at the root of the 35° swept wing used in 
the present investigation is shown in figure 12 . The effects of taper 
and aspect ratio for this wing are believed to he small and were neg- 
lected. The slope of the airfoil surface was determined graphically. 
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Body Modification 


The "body is replaced by a vortex cylinder, that is, a cylinder upon 
which ring vortices are distributed so as to induce an axial -velocity 
ratio Um/Vo of sufficient strength to cancel the velocity ratios 
Ub/V 0 and u c /V Q . The vortex cylinder is assumed to be of constant 
radius, rj, equal to the body radius at the midpoint of the reference 
chord (at the wing -body junction) and of length, c j, equal to the length 
of the reference chord. (See fig. 3*) 

Six standard distributions of vortex rings, suitable for the design 
of wing -body intersection lines, and the corresponding induced velocities, 
are given in reference 7» The distribution functions are: 

Even functions 

7 s - = </l - (21 - l) 2 

/o 2 

7 s - = V'l - (21 - l) 2 
°6 

7 s - = l - (24 - l) 2 
/o 8 

Odd functions 

T 2 - = - (21 - 1) J 1 - (24 - l) 2 

°3 


—■ = - ( 24 - 1 ) yiT(24 - 1) 2 

/q 7 


7 s - = - (24 - 1) 1 - (24 - l) 2 

7 0 9 L 


Here, a coordinate system is used in which 4 is made dimensionless 
with respect to the reference chord, positive when directed downstream 
and with origin at the leading edge of the reference chord. The func- 
tions 7y/h 0 are plotted in figure 13 * 
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The distribution of ring-vortex strength along the cylinder is 
determined by equating the induced axial velocity ratio to the value of 

^ at a number of control points along the reference chord. 

The corresponding radial component of induced velocity is then used to 
determine the body- shape modification Ar. 

The axial velocity induced by the vortex cylinder can be repre- 
sented by 




(A7) 


and the radial increment Ar by 



(AB) 


where 7 v /7 0v , u v */v o , and Ar v */(cj/2) are functions tabulated in 

tables III, IV, and VI, respectively, of reference 7 for various values 
of the vortex- cylinder length-to-diameter ratio. (Values of x in the 
tables of reference 7 correspond to values of 2£ - 1 in the notation 
used here.) 


Since there are six vortex distribution functions to be used, the 
velocity ratio u m /v o must be evaluated at six chordwise control points. 
For the work of this report, control points at i = 0.1, 0.3, 0.5, 0.6, 
0.7, and 0.9 were chosen. The tables in reference 7; which are for 


M = 0, 


may be used directly if the velocity ratio 


( 


u b + u c \ 

Vo / 


is reduced 


"by 3 des and if a vortex cylinder of reduced radius 3 des rj is used. 

(See ref. 14.) Equation (A7) may then, with the help of tables III and IV 
in reference 7, be applied at each control point and the resulting system 
of six equations can be solved simultaneously for the six unknown coef- 
ficients C v . When the coefficients Cv are known, the modification Ar 
can then be calculated by using equation (A8), for which tabulated values 
of Ar v /(c j/2) are to be found in table VI of reference 7* 
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As an illustration, the radial modifications for the body having a 
fineness ratio of 12 (wing in the forward position) are shown in fig- 
ure l4(a) . It is readily seen that the radial modification (calculated 
by summing the odd and even functions) rapidly approaches zero upstream 
and downstream of the reference chord. For application in the present 
method of body contouring, the value at I = -0.2 was subtracted in 
order to start the modification at that point, and an arbitrary fairing 
was used to terminate the modification at i = 1.2. (See fig. 14(b).) 


Additional Modification 


In the preceding development, the vortex intensity was constant 
around the periphery of each cylindrical cross section, so the calcu- 
lated changes in body radius resulted in annual bumps, invariant with 
position around the circle. The major portion of the "interference 
velocity" u m , however, is attributable to the wing and attenuates 
rapidly in the vertical, direction. An approximation for the attenuation 
at the top of the body was obtained by considering the variation with 
vertical distance of the maximum perturbation velocity induced by a two- 
dimensional, unswept, biconvex airfoil. According to linearized theory, 
this variation is 



With this expression as a guide, the radial modification at the top of 
the body was reduced by the amount 


Ar z 

Ar 


= 1 - 


2Pdes r j 


tan -1 


20 


'des r i 


(A10) 


and the contoured body shape was made elliptical in cross section with 
the full radial increment at the side, so that Ary = Ar. Details of 
the contoured bodies are presented in figure 3. 
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A-16454 

(a) General view of swept wing mounted on a modified small body having a 

fineness ratio of 9* 



A-16451 

(b) General view of swept wing mounted on a modified large body having a 

fineness ratio of 9* 

Figure 1.- Photographs of the models. 
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(c) Close-up of intersection of swept wing with modified large body hav- 
ing a fineness ratio of 9* 



(d) Close-up of intersection of swept wing with basic large body having 

a fineness ratio of 9* 


Figure 1.- Concluded. 
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(a) Body having a fineness ratio of !2. (See fig. 2(a).) 
Figure 3 -Contouring details . 
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(b) Small body having a fineness ratio of 9. (See fig. 2(d).) 


Figure 3 -Continued. 
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(c) Large body having a fineness ratio of 9. ( See fig. 2 (cl.) 
Figure 3 -Concluded 
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Mach number, M 

Figure 4.- Variation of Reynolds number with Mach number. 
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o Basic body 
□ Cylindrical body 
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(b) Small body with a fineness ratio of 9. 
Figure 5 - Continued. 
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U 




Figure 6 -Variation of drag coefficient at zero lift with Mach number 
for the swept wing in combination with various bodies. 
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.5 .6 .7 .8 .9 i.O 

Mach number, M 

(c) Large body with a fineness ratio of 9. 

Figure 7- Variation of drag coefficient at 0.3 lift coefficient with 
Mach number for the swept wing in combination with various 
bodies. 
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( a) Body with a fineness ratio of !2 



( b) Small body with a fineness ratio of 9. 



Mach number, M 

fcj Large body with a fineness ratio of 9. 

Figure 8 -Variation of lift -curve slope with Mach number for the swept 
wing in combination with various bodies. 
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(a) Body with a fineness ratio of 12 
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Mach number, M 

(c) Large body with a fineness ratio of 9. 



Figure 9 .-Variation of stability parameter with Mach number 

l C l =0 

for the swept wing in combination with various bodies. 
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o Measured 
— Calculated 



(a) Body having a fineness ratio of !2 (wing forward). 



(b) Body having a fineness ratio of !2 (wing aft). 


-.12 
-.08 
-.04 

0 

!2 16 20 24 28 32 

(c) Small body having a fineness ratio of 9. 




(d) Large body having a fineness ratio of 9. 


Figure //.- Pressure distribution from body- alone tests compared 
with theoretical pressure distributions calculated by equation 
(A3) for a Mach number of 0.87. 
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Figure 12.- Distortion velocity at the root of a 35° swept back, un tapered 
wing having 64 Z A0I5 sections normal to the leading edge (Mach 
number of 0.87). 
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Figure i 3 -Vortex distributions. 
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(a) Radial modification for the odd and even vortex- 

distribution functions. 



e 


(b) Radial modification adjusted to have zero value at 

£= - 0 . 2 . 

Figure 14- Radial modification as computed by the ring- 
vortex method for the body having a fineness ratio 
of !2 (wing in forward position). 
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